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We consider the electromagnetic waves propagating in the system of coupled waveguides. One of 
the system components is a standard waveguide fabricated from nonlinear medium having positive 
refraction and another component is a waveguide produced from an artificial material having nega- 
tive refraction. The metamaterial constituting the second waveguide has linear characteristics and 
0^ ' a wave propagating in the waveguide of this type propagates in the direction opposite to direction 

' of energy flux. It is found that the coupled nonlinear solitary waves propagating both in the same 

' direction are exist in this oppositely-directed coupler due to linear coupling between nonlinear pos- 

itive refractive waveguide and linear negative refractive waveguide. The corresponding analytical 
^ ' solution is found and it is used for numerical simulation to illustrate that the results of the solitary 

wave collisions are sensible to the relative velocity of the colliding solitary waves. 
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I. INTRODUCTION 



The waveguide structure fabricated from the two closely placed waveguides is of common use in the fiber and integrated 
^ optics. Coupling between the waveguides is due to tunnel penetration of light from one waveguide into another 
__i waveguide [l], Q • This coupler preserves direction of light propagation, and for this reason it is named a directed 
coupler. It was found [1, that the steady state pair of electromagnetic pulses can exist in the extended directed 
~^ ' coupler or twin-core fibers § • Sometimes it termes soliton. There is a large body of publications devoted to 
^ , investigation of the soliton generation and its propagation, see for example @, 



' Recent prog ress in nanofabrication has lead to the design of new materials with highly unusual optical properties 



[Til. [T^. ITsf. 1 14l. ITsl ITsl. [TqI. [2^ . of whichnegative refraction is an example. Negative refraction occurs in media in which 



§ . the wa.e vec.» o, the e,ect„„,ag„etic wave is a.,«p»al,e, to the Poy„t.„g vecto. M M the partiCa. eaee 

■ ' ' when the real parts of the dielectric permittivity and magnetic permeability in the medium simultaneously take on 
negative values in some frequency range, the property of negative refraction will appear. The existence of such media 
Qs^ was demonstrated experimentally first in the microwave and then in the near-infrared ranges [H, [H, H) [Hj [S [13 ■ 
I Negative refraction can be realized for the media with spatial dispersion 20] which is large enough and in the photonic 
crystals [24, 25, 26]. Recently bulk metamaterials that show negative refraction for all incident angles in the visible 
region were presented 0, Hsj . 

' The unusual properties of negative refractive (NR) index materials reveal themselves most prominently when the 



J-j ' refractive index of the same medium can be positive in one spectral region and negative in another |23l . |29| . l3d . |3l| . 
. . ■ New wave propagation phenomena can also be expected when a wave passes through, or is localized near, an interface 



between such a material and a conventional dielectric, i.e., positive refractive (PR) index material [32, l33|, |3J, |35| • We 
can refer to these cases as negative- positive refraction ones. The intriguing example of the negative-positive refraction 
medium is the coupler, where one of the waveguides is fabricated from a material with a negative refractive index. 
This device acts as a (distributed) mirror. The radiation entering one waveguide leaves the device through the other 
waveguide at the same end but in the opposite direction. For this reason, this device can be called a oppositely- 
directed coupler (ODC). The ODC is known to support the propagation of linear waves with a gap in their spectrum 
(forbidden zone). It has been shown recently 37] that due to this gap the nonlinear oppositely directed coupler is 
bistable. Bistability results from the multi- valued dependence of the transmission coefficient on the input-wave power. 
It is noteworthy that this effect has no analogy in conventional directed couplers consisting of uniform waveguides 
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without a mirror-based feedback mechanism. It was found ^3§\ that if the waveguide nonhnear optical properties are 
characterized by the third order susceptibihty, a coupled pair of the steady state pulses (each one is localized in its 
own waveguide) can exist in ODC. 

The fabrication of the transparent nonlinear NR waveguide is a complicated and still unresolved problem. The 
linear NR waveguide could be more suitable for realization. However there is question is it possible a quasisolitonic 
regime for the waves in such a antisymmetric coupler consisting of a nonlinear conventional dielectric waveguide 
coupled with a linear NR waveguide. In this paper we consider the extended asymmetric oppositely-directed coupler 
(AODC). We found the steady state solutions of the equations which describe the waves in the AODC. The numerical 
simulation of these steady state pulses formation and the collisions between the quasisolitons is performed. We found 
that the pulses are very robust against perturbations. It allows a definite conclusion that the steady state pulses are 
a new kind of gap solitons. It should be noted that this solitary wave exists in the medium without grating. The 
influence of linear losses in the NR waveguide on the quasisoliton formation discussed shortly. 



II. THE MODEL FORMULATION 



We consider the waveguide structure in which one of the channels is made from nonlinear material features by the 
positive refraction, and the NR channel is composed of linear material. The system of the equations based on connected 
wave approximation takes the form 

dqi , i dqi ^ .^Af3z , 



i-^ + + -^12926 

az Vgi ot 



27rt^o / mi(^o) ^2, ,(3) 
c y ei(wo) 



^oxrAlgilV^o, (1) 



dq2 i dq2 _,ai3z n 

OZ Vg2 ot 

The coefficient xfff effective non-linear susceptibility of the first channel. K12 and K21 are the coupling 

constants. Linear properties of the PR channel defined by dielectric permittivity ei and magnetic permeability /ii. 
These values are assumed be real, that corresponds to the lossless materials. Slowly varying envelope of the electric 
field Ej{z, t) in j-th channel is written as Ej{z, t) = AQqj{z, t), Vgj is the group- velocity of the wave in j'-th channel, 
(j — 1, 2). Here we assume that the group- velocity dispersion can be neglected. In the following we also assume that 
the synchronism condition is satisfied: A/3 = 0. It is convenient to use new normalized variable 

Oi = VK^ilie-''^^^ Q2 = /?^g2e'^''^ 
C = z/L„ T = t^,\t- z/Vo), L, = {Ki2K2i)-^'^, 

ta = Lc{Vgl + Vg2)/2VglVg2, Vq^ = [v gl - U<,2 ) / 2 1 W (,2 • 

In the normalized variables the ultimate system of equations ([T]) for the AODC reads as 



{■k + ^)Qi + Q2+r\Qi?Qi = 0, 



(2) 



The parameter of nonlinearity is 



27ru;o / a^i(^o) .2 (3) 



CK21VK12K21 V ei(wo) 
Let us consider the linear wave limit of the AODS equations. From ^ one can get 



If we take the Fourier transformation 



9C 

i ~ Q2 — Qi = 0. 



(3) 



1,2 = / Qi,2e 



4772 ' 
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the equations ([3]) result in the following linear system of equation 

(i/ - k)Qi + Qa = 0, 
{iy + K)Q2 + Qi =0. 

This system of equations has the nonzero solution only if the corresponding determinant 

detf'^T^ ! 
is equal to zero. That leads to the dispersion relation 



=±v/1 + k2. (4) 

Thus the spectrum of the linear waves has the gap Avg — 2. This gap is characteristic feature for a distributed mirror 
[3^. Hence, the AODC in linear wave limit acts as a mirror. 

In following consideration it is suitable to take the real variables form of equations By using the real variables 
Qi,2 — 01,26*'^^'^ one obtains 

_ + _j,,^,,sm$, (^_-_ja,.a,smci>, 
d d \ I a2 ^ , 2 



7^ + 7r '^1 = — cos$ + raf, (5) 



- TJ- 02 = COS$, 



where ^ — cjji — (f)2- From the amplitude equation it follows that 

d , n n^ d 



[a 



Hence, 



d_ 



+ 00 

2\ 1+^ 



{4~al)dr^ {al + ai)\ZZ- (6) 



In the case of solitary waves, for which the electromagnetic fields vanish at infinity are considered, the right part of 
^ is equal to zero. It leads to integral of motion 



+00 

,2 2 



a2 — fli) dr — const. (7) 

It is the modified Manley-Rowe relation. As usually, in the case of the quadratic nonlinear PR media Manley-Rowe 
relation looks like 

+00 



J + af) dr — const. 



The difference of the expression ([T]) from conventional Manley-Rowe relation is explained by the fact that the energy 
flux for the waves in one waveguide opposite in direction to energy flux of other waveguide, while their wave vectors 
are approximately the same. This flux pattern is an inherent feature of the negative-positive refraction media. 
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III. ANALYTICAL SOLUTION — THE STEADY-STATE PULSE 

To consider the solitary steady state waves in AODC by similar way as in jsl] we have start from the equations ([5]). 
Suppose that solutions of these equations are depend only on single variable 

V 



with free parameter (3. Suppose ui = yT+TJoi and U2 = ^/l — (3a2. The system of the equations ^ takes the 
following form 

d d 

— Ui=M2sin$, — U2 = uisin$, (8) 
arj orj 

= —cos $ + 9^2, (9) 
orj ui 

^ i2 = -— cos$, (10) 



where 



drj U2 



e 



'1-/3 



1 + I3\] l + f3 

We can also write an equation for the phase difference 



(^!^ + !^ )cos$ + eH^ (11) 



dr] \U2 ui 

The phase equations can be used to get one equation for the phase difference. Finally, the total system of equations 
reads 

d 

—Ml = U2sin$, (12) 
or] 

d 

—U2 = Misin$, (13) 
orj 

— $ ^ _L + _i cos$ + eu?. 14 
or] \U2 uij 

We are looking for a solution in a form of the solitary wave, it corresponds with the following boundary condition 

ai,2 ^0 at r] ±oo. 

From the equation (jl2|) and (fT3|) it follows u\ ~ u\, or ui — eu2, where e — ±1. Hence, the system of equations 
(fT2|) - (|14p is reduced to following pare of equations 

d 

— Ui = euisin<i>, (15) 
or] 

d 

— $ = 2ecos<I> + eu?. (16) 
orj 

Multiplying the last equation by sin$ and taking into account the equation (|15p . we get the second integral of 
motion 

ul (^cos<i> + ^ul^ =C2. (17) 

Due to the boundary condition, the value of this integral is equal to zero. As we are looking for a non zero solution 
of (fT5)) and (I16p . one can write 

cos$ + ^mJ = 0. (18) 
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Substitution of ([18)) into ((15)) leads to 

{du,/drif ^ui{l- {Q/Afu\) . 

Choosing the variable ui — w^^^"^ reduces this expression into equation 

{dw/df]f =4(u;2_(e/4)2)^ 

which has the following solution 

■u;(?7) = (6/4) cosh 2(?7- rya). (19) 

Thus, the solution of ()T5)) and ((TC)) is 

2, ^ _ 2,„^ 4/e 



cosh 2(77 - 772) ' 
By using the phase's equations and (fTO)) one can write 

—01 = -eu2(?7) 



9?? 4 cosh 2(77 — 772) ' 

9 1 1 

— (/)2 = -tQuHv) = — X- 

Of] 4 cosh 2(77 — 772) 

that yields 

01(77) = 0i(-oo) + 3arctan(exp2(7; - 772)) , (20) 

Mv) = 02(-oo) + arctan (exp2(77 - 772)) . (21) 

The phases 0i_2(— 00) should be chosen in such a way that (llSp is satisfied. As the solitary wave's amplitude tends 
to zero ui,2(oo) 0, the relation p8|) at infinity reduces to 

cos $(—00) — cos(0i(— cx)) — 02(— 00)) = 0. 

In the numerical simulation we choose 0i(~oo) = and 02(~oo) = — 7r/2. 
The amplitudes ai,2 are defined by the following expressions 

4 

"'^"^^ " 6(1 + /3) cosh2(77-7y2)' ^^^^ 

"'^'^^ " 6(1-/3) cosh 2(77 -772)- ^^^^ 

The expressions (|22)) and (|23)) describe the steady state solitary wave propagating in AODC under consideration. 
This wave looks like a gap soliton propagating in the nonlinear Bragg grating, however there are no any periodic 
structures. Gap in the linear wave spectrum (jj)) is due to the flux pattern in this negative-positive refraction media. 

Some remarks are worthy to be made about the free parameter (3. The negative value of the parameter /3 corresponds 
to the solitary wave propagating in the direction of the axis C- The solitary wave characterized by positive value of the 
parameter f3 propagates in the opposite direction. The large amplitudes (more powerful solitary waves) correspond 
to large positive values of the parameter /3. For the negative values of the parameter /3 the quasisolitons with smaller 
amplitudes have smaller values of the parameter (3 however the absolute value of the velocity determined by the 
parameter (3 is larger for the less powerful solitary waves. We will refer to the Fig[T] later in discussion the influence 
of the linear losses at the quasisoliton. 



IV. NUMERICAL SIMULATION 



It is common knowledge that the completely integrable evolution nonlinear equations have a special solutions describing 
elastic interaction of solitary waves [40|. These waves are named solitons. It is our opinion that the system of equations 
(12) does not belong to the class of completely integrable equations. Hence the solution of these equations does not 
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FIG. 1: (Color online)The plot for the solitary wave amplitude ai,2 dependence on parameter /3. Red solid curve corresponds to the 
amplitude of the pulse in the nonlinear PR channel, dashed green curve is for the NR channel. 




10 20 30 40 50 60 10 20 30 40 50 60 

T=(t-z/Vo)/to x=(t-z/Vo)/to 



FIG. 2: Crossing collision between two solitary waves with (i = —0.9 and /3 = —0.5. Left panel is for the solitary wave in the 
PR waveguide, right panel is for the solitary wave in the NR waveguide. 

represent true soliton. We will denote them as quasisolitons. To investigate interaction between the steady state 
solitary waves ()22|) and (|23|) the numerical simulation was pursued. 

Numerical solution of the system of equations ^ is performed using the scheme of the finite differences. The 
conservation of the first integral of motion ([7]) is used to control the computational error. In the consequent an 
integration step over the evolution variable r is set to hr = 0.0005. The nonlinearity is set to r — 1. (The smaller 
nonlinearities correspond to the quasisolitons with the same profiles but the smaller amplitudes (defined by the 
parameter O which is linearly proportional to the nonlinearity coefficient.) 

It was found that the results of the collision depend on the relative velocity of the pulses. Fig[2] illustrates collision 
between quasisolitons characterized hy (3 = —0.9 and (3 = —0.5. The relative interaction velocity is 0.4. The 
quasisolitons drop some radiation after interaction and weakly radiating wave could be noticed in the left panel of 
the picture as a result of collision. Velocities and the amplitudes of the resulting quasisolitons are slightly different 
from the parameters of the quasisolitons before the collision due to loss of some part of energy. 

Decrease the difference between velocities of colliding pulses results in long-range interaction between pulses. The 
interaction distance depends on the relation of the velocities. In the FigO is shown interaction of two quasisolitons 
characterized hy (3 — —0.7 and (3 = —0.5. The relative interaction velocity is 0.2. The quasisolitons exchange energy 
at the distance and form a transitory bound state. 

It was mentioned before that the quasisolitions with the same absolute values of velocity have different amplitudes. 
To investigate the robustness of the quasisoliton depending on its energy, at the Figl3]we present results of modeling 
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FIG. 3: Repulsive collision between two solitary waves with /3 — —0.7 and /3 = —0.5. Left panel is for the solitary wave in the 
PR waveguide, right panel is for the solitary wave in the NR waveguide. 
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FIG. 4: Collision between two solitary waves with /3 = 0.7 and f5 = —0.7. Left panel is for the solitary wave in the PR 
waveguide, right panel is for the solitary wave in the NR waveguide. 



the collision between two steady-state solitary waves with the same absolute values of the velocity (3 — 0.7 and 
(3 — —0.7. A more energetic quasisoliton with /3 = 0.7 remains unchanged after coUision and the less energetic 
quasisoliton with [3 = —0.7 loose some radiation and changes its trajectory. 

Thus, the coUision of two steady state pulses with different velocities has shown significant robustness. SmaU 
amplitude radiation appearing after coUision attests that the quasisohton eventuaUy be disappeared. 

For the current state of fabrication technology the losses in the real NR materials are considerable. High value of 
losses renders steady state pulse propagation impossible. Nevertheless, the role of small losses would be considered. 
To check the influence of linear losses in the NR channel on the solitary wave formation and propagation, the model 
equations ([2]) were modified by including the additional term in the equation for a NR channel 

« (^^ + ^) 0i + 02+r|Qi|2Qi =0 (24) 

^(|-^-7)Q.-Qt =0. (25) 

The FiglS] illustrates the evolution of same pulses as in the FigH] placed at the boundaries however the linear losses 
in the NR channel are taken into consideration 7 = 0.05. Comparing these two pictures one may conclude that 
even small losses affect considerably the propagation properties of the solitary waves. Propagating in the dissipative 
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FIG. 5: Two solitary waves with /3 = 0.7 and f) = —0.7 in a system with the lossy NR channel 7 = 0.05. Left panel is for the 
solitary wave in the PR waveguide, right panel is for the solitary wave in the NR waveguide. 

medium these pulses loose energy and the pulse with positive /3 slows down to sustain the steady state regime. The 
pulse with negative (3 slightly accelerates to support the quasi-steady regime of propagation (notice a Fig[T]). 

V. CONCLUSION 

We considered the nonlinear solitary waves propagating in the nonlinear opposite-directed coupler. One of its compo- 
nents is a nonlinear waveguide made from the material with positive index of refraction. Another channel is fabricated 
from linear dielectric material with negative refraction index. For the system of coupled waves the analytical solution 
for the electromagnetic waves in the form of steady state solitary waves is found. It is interesting to notice that the 
wave in the nonlinear channel affects the wave propagation in the neighboring linear channel. This influence results 
in coupling and a steady propagation of the solitary waves in both waveguides, the wave in the nonlinear waveguide 
draws the wave in the linear NR waveguide. 

It was shown that the result of interaction of the quasisolitons depends on their velocities. Provided that relative 
velocity of the two colliding quasisolitons is large, the pulses collide almost like the solitions. A weak radiation or a 
weak pulse emerges as a result of collision. At small relative velocities of the quasisolitons the interaction between 
two colliding solitary waves could results in a strong energy exchange between colliding solitary waves and formation 
of the temporarily coupled state of the interacting quasisolitons. This phenomena is known for the case of nonlinear 
directed coupling. In the case under consideration the formation of the coupled state demands for an additional 
investigation. It could be possible that a long-living coupled state of two solitary waves is exists there. 

The influence of the linear losses in the NR waveguide on the existence of the solitary waves was studied. The wave 
loses its energy as it propagates through the dissipative medium. In order to compensate the losses the entering pulse 
in a form of the solitary wave characterized by positive value of parameter (3 (its propagation direction is opposite to 
the coordinate axis C) permanently transforms into another solitary wave, less energetic, having smaller value of the 
parameter (3. It results in slowing down the wave in the lossy medium. After that the wave stops in the waveguide and 
changes the direction of its propagation. Oppositely the entering pulse in a form of the solitary wave characterized 
by negative parameter /3 (its propagation direction is coincides to the coordinate axis C) tends to increase its velocity 
(i.e. to decrease parameter (3 which is negative) to sustain the steady-state regime of propagation in presence of the 
losses. 
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